Onset of ferromagnetism in low-doped Gai_xMnxAs 
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We develop a quantitatively predictive theory for impurity-band ferromagnetism in the low-doping 
regime of Gai_xMnxAs and compare with experimental measurements of a series of samples whose 
compositions span the transition from paramagnetic insulating to ferromagnetic conducting behav- 
ior. The theoretical Curie temperatures depend sensitively on the local fluctuations in the Mn-hole 
binding energy, which originates from disorder in the Mn distribution as well as the presence of 
As antisite defects. The experimentally-determined hopping energy at the Curie temperature is 
roughly constant over a series of samples whose conductivities vary more than 10* and whose hole 
concentrations vary more than lO'^. Thus in this regime the hopping energy is an excellent predictor 
of the Curie temperature for a sample, in agreement with the theory. 
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The III-V ferromagnetic semiconductor Gai_xMnxAs 
has been the focus of intense interest due to its potential 
incorporation in proof-of-concept spintronic devices [TJ |5J 
[31 H] . Most studies have focused on relatively high Mn 
concentrations (x ~ 0.03 - 0.08) [3 H H [8. for which the 
Curie temperature can reach values of well over 100 K. 
Theoretical treatments of the Curie temperature in the 
low doping regime [21 UHl [HI [Hj have focused on the 
Mn dopants in the sample, including their concentra- 
tion, distribution, compensation, and whether they are 
substitutional or interstitial. Those that treated variable- 
range hopping (VRH) [TT| considered the Efros-Shklovskii 
regime [13], where the hopping energy arises from the 
Coulomb interaction, and found an even stronger depen- 
dence of the Curie temperature on hole doping than in 
the higher-doped regime. Recent experimental measure- 
ments on low-doped Gai_xMnxAs, however, show the 
dependence of the Curie temperature on hole doping to 
be negligible, or even nonexistent [TJJ , such that samples 
with the same measured hole density and Mn density may 
differ in Curie temperature by a factor of 2. This sug- 
gests the proper VRH regime is Mott VRH [HI [Ml [TT] . 
in which hopping energies are independent of the hole 
density. 

Here we present a quantitative theory of the Curie 
temperature of Gai_xMnxAs in the Mott VRH regime 
which is built on a microscopic model for the effect of 
disorder on spin-spin interactions. We also present mea- 
surements of the Mn concentration, hole concentration, 
conductivity, hopping energy, and Curie temperatures in 
a series of samples prepared via a combinatorial approach 
that exploits a spatial gradient in the Ga:As flux ratio 
on non-rotated substrates during MBE growth. We use 
our theory to estimate the hopping energy from the den- 
sity of substitutional Mn (ionized or neutral), and find 
that it is comparable to the ferromagnetic interaction en- 



ergy between two spins. We also present self-consistent 
calculations of the Curie temperature in the Mott VRH 
regime, and we find the Curie temperature to be con- 
trolled primarily by the hopping energy, in agreement 
with the experimental findings. 

We begin by describing our theoretical model and re- 
sults. The typical hopping distance and hopping energy 
play an important role in Mott VRH, for sites that are 
spatially close typically have energy differences too large 
to permit rapid hopping, whereas those that are spatially 
far away have small hopping matrix elements. The hop- 
ping energy W is the typical energy difference between 
two sites separated by the most likely hopping distance 
£ of a carrier, and can be directly extracted from the 
temperature dependence of the conductivity, which fol- 
lows a T^l^ dependence. As noted in Ref. [TT], similar 
quantities influence ferromagnetic interactions between 
Mn dopants. In our model, however, W and £ are deter- 
mined by the availability of states, rather than interac- 
tions between holes. Shown in Figure [T](a) is a schematic 
double exchange model for two interacting Mn with par- 
allel spins, including a site energy difference Ai?. As in 
typical two-site models of ferromagnetism jl 81 , for paral- 
lel spins the hole bound state associated with Mn #1 
hybridizes with that of Mn ^^2, and one molecular state 
has lower energy than either isolated Mn state. If a sin- 
gle hole is present per pair then parallel alignment of the 
Mn spins is preferred to antiparallel alignment (in which 
the states cannot hybridize). Bound hole states at dif- 
ferent energies A_B hybridize more weakly and lead to a 
reduced ferromagnetic coupling as shown in Fig.[ljb). 

We calculate the interactions between two Mn in the 
presence of disorder, H = J Si -82, in this dilute car- 
rier regime. Shown in Fig. [l] are calculated values for 
J{AE) / J(0) as a function of AE obtained from a multi- 
band tight-binding theory [12] that has successfully pre- 
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FIG. 1: (color online) (a) Model taking into account an en- 
ergy difference /S.E between the two Mn hole states, (b) Fer- 
romagnetic interaction energy for Mn-Mn pairs separated by 
17A. Line between the two Mn is along the [221] direction, 
with sample magnetization along [001] (red) or [100] (blue), 
or along the [100] direction, with sample magnetization along 
[001] (orange) or [100] (green), (c) Tc as a function of To for 
several averag ed values of 27(0)5"^ for 17A pairs. 

dieted the interaction energies of Mn-Mn pairs measured 
by STM[20]. J(0) depends on the specific Mn-Mn pair, 
yet the dependence of J{AE) on AE is reproduced by 
an analytic expression from the model in Fig. [ija), 

J(A£;)/J(0) = {{AE/AJ{0)S^f + iy^^-\AE\/AJ{0)S^. 

(1) 

Plotting the results for Mn-Mn pairs separated by 17A 
in scaled units collapses all the curves onto the universal 
shape shown in Fig. nib). At a distance of 17A the calcu- 
lated interaction energy 2J{0)S^ varies from 8.5 meV to 
36 meV, depending on the pair geometry. For other dis- 
tances the average interaction energy can be estimated 
using the approximate decay length of the acceptor wave 
function, = 13A.|2T] 

To evaluate the effect of hopping transport on the fer- 
romagnetic interaction we assume the most important 
pairs for carrier-mediated ferromagnetism are those pro- 
viding the dominant contribution to the conductivity, 
and use the hopping energy 

WiT) = [T{5/A)/2]kBT^^^T^/^ (2) 

for AE. From the results of Fig. [T] once W becomes of 
the order of 2J{0)S^ then the ferromagnetic interaction 
will be significantly reduced. We follow Ref. [TUto evalu- 
ate i and Tq. Calculating the hybridization between Mn 
sites [m [2^ separated by the typical inter-Mn distance 
for X = 0.015, I5A, yields an estimate of 200 meV for the 
bandwidth, similar to that found in Ref. [TT] There are 
3 states per Mn, and thus /cbTq ^ 700 meV. The Curie 



temperature can be found by solving 

iiTc) = a-ir(5/4)(ro/Te)i/V2, (3) 

Tc = S{S+l)MW{Tc))/3kB. (4) 

simultaneously with Eqs. ([l])-([2|. T^s are shown in 
Fig.[l];c) for several values of 2J{0)S^ for 17A pau's. 

We now compare the results of this model with the 
magnetic and electronic properties of a series of sam- 
ples within which the electronic properties were continu- 
ously varied as described in detail in Ref. [14]. We grew 
100 nm thick Gai_xMnxAs epilayers on top of a 300 nm 
GaAs buffer on semi-insulating (001) GaAs wafers for x 
= 0.0075, 0.01, 0.0125, and 0.015, where the quoted Mn 
concentration values are calibrated from reflection high- 
energy electron diffraction (RHEED) measurements and 
secondary ion mass spectroscopy (SIMS). Due to the ar- 
senic flux gradient [13] , the carrier density varied across 
each wafer by at least two orders of magnitude from com- 
pensation by As antisites. We cleaved three long strips 
side-by-side (~ 4 mm in width) from the center portion 
of each wafer along the As-gradient direction. Each strip 
was then cleaved into 15 sample pieces (~ 3 mm x 4 mm) 
for transport, magnetization, and SIMS measurements of 
samples with varying As:Ga stoichiometry. 

We measured magnetization with a commercial super- 
conducting quantum interference device (SQUID) mag- 
netometer (Quantum Design MPMS), whereas resistivity 
and Hall effect measurements were performed in com- 
mercial cryostats (Quantum Design PPMS) using four- 
probe techniques with external electronics to measure re- 
sistances up to > ICri. We determined the existence 
of ferromagnetism, and Tc, from temperature-dependent 
magnetization data, taken on warming in a 50 Oe field af- 
ter initial cooling from room temperature to 2 K in a 1 T 
magnetic field. Measurements of M{H) hysteresis loops 
confirmed the existence of ferromagnetism. We used 
annealed indium contacts on lithographically patterned 
Gai_xMnxAs hall bars for transport measurements, and 
the contacts were confirmed to be ohmic in the source 
current range of 10 - 1000 uA and T = 2 - 300K. 

Figure [2] shows the temperature dependent resistivity 
p{T) and the magnetization M{T) along the [100] crys- 
talline direction of Gai^xMuxAs epilayers synthesized us- 
ing the conventional growth technique (rotating the sub- 
strate for uniformity). The peak observed near Tc in 
the p{T) data of Fig. [2]ja) is typical behavior for metal- 
lic high-Tc Gai_xMnxAs.[Tl |2] For smaller values of x, 
near the onset of ferromagnetism, insulating behavior 
(dp/dT < 0) is observed at all measured temperatures, 
as shown in Fig. [2f^b), and there is effectively no feature 
near Tc- These data clearly demonstrate a ferromagnetic 
insulating state in Gai_xMnxAs. 

We interpret the temperature dependence of the resis- 
tivity through the VRH expression p = pa exp{TQ/Ty/^ , 
although we caution that this analysis can only be ap- 
plied to our data over a limited temperature range. 
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FIG. 2: (color online) Temperature dependence of the resis- 
tivity (open symbols) and the magnetization (closed symbols) 
along the [100] crystalline direction of as-grown (Ga,Mn)As 
epilayers grown in rotated mode for x — 0.03 and 0.015. 



W{Tc) is then calculated from Tq via Eq. Figure[3]^a) 
shows Inp(r) as a function of T~^^^ between 80 K and 
25 K (i.e. above T^) for all of the pieces from the x 
= 0.015 non-rotated sample. The solid lines are linear 
fits performed in this temperature range, and the single 
dashed line in the center of the graph separates the ferro- 
magnetic (FM) and non-ferromagnetic (non-FM) regime 
(as determined from magnetization measurements). Sim- 
ilar trends and quality of fit were observed in the trans- 
port data sets for the other Mn concentrations studied. 
Figure [sjb) shows the full temperature range of 2 — 300 K 
for all FM samples for the x = 0.015 series. A noticeable 
deflection appears near Tc, suggesting that the magnetic 
and transport properties are still connected even in this 
low doping regime (although a similar deflection has been 
reported in a non-magnetic doped GaAs system) |22J . 

We plot the Tc of all of our non-rotated samples in Fig- 
ure |4] as a function of carrier concentration and conduc- 
tivity at T = 25 K, and hopping energy at the Curie tem- 
perature. The carrier concentration is derived from Hall 
effect measurements between —2 T and 2 T taken at 300 
K due to the high resistance at low temperatures. Every 
data point in Fig. [4] represents an individual cleaved piece 
from the wafer strips as described before. Although not 
explicitly shown, there are no conductivity data points 
for W{T = 25K) > 21 meV since the sample resistance 
at 25 K is too large to measure with our apparatus. De- 
spite the exceptionally broad distribution of conductivi- 
ties (over four orders of magnitude) and hole concentra- 
tions (over two) , all of the ferromagnetic materials have 
W{Tc) between 3 and 10 meV. 

The value of £ extracted for the ferromagnetic samples 
is 22A (as i depends on Tc^^, there is little variation 
in i for this range of samples). The estimated value of 
27^(0)5'^ for this distance is then 8 meV, which is within 
a factor of two of the measured hopping energy at the 
Curie temperature for all the samples, in agreement with 
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FIG. 3: (color online) Resistivity from 25 K to 80 K for x 
— 0.015 non-rotated samples. The straight lines of the same 
color are the corresponding fits performed in this temperature 
range as discussed in the text, (b) The full temperature range 
T = 2 - 300 K for all FM samples from the x = 0.015 series, 
where a noticeable deflection appears near Tc and separates 
the curve into two segments. 



the theory. The bandwidth, which influences Tq, does 
depend on the Mn concentration, but for a change of Mn 
concentration of a factor of 2 both W and £ will only 
change by 2^/^, which for the measurements made here 
has a negligible effect on Tc- The absence of ferromag- 
netic samples of < lOK is consistent with the solutions 
shown in Fig.[ljc), where a discontinuous drop from a fl- 
nite Tc to non-ferromagnetic material is seen. 

The high carrier density limit (~ 10^^ cm~'^) of these 
low Mn-doped samples just reaches the lower limit of 
metallic Gai_xMnxAs considered in typical theoretical 
calculations. Compared to other Gai_xMnxAs samples 
with a similar Mn-doping level, the resistivity of our fer- 
romagnetic samples is several orders of magnitude lower, 
which is possibly due to an optimized control of the ar- 
senic and the reduction of As-antisites. This reduced 
disorder will lead to a smaller bandwidth, and thus a 
smaller Tq, which (see Fig. [ijc)) permits the materials 
to be ferromagnetic. Our data show that the carrier 
concentration has a weak effect on the onset of ferro- 
magnetism in this limit. Instead, the strong correlation 
of the cutoff in ferromagnetism with critical values of the 
electric conductivity and the hopping energy suggest that 
the interplay between disorder and Mn-Mn interactions 
found in Eqs. ([l])-([4]), and shown in Fig. [ijc), dominate 
the transition. The data in Fig. [4] encompassing sam- 
ples over a broad range of As:Ga stoichiometry, indicate 
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magneto-optoelectronic applications. The theory pre- 
sented here might also be extended to higher hole den- 
sities. Screening by holes of disorder potentials would 
lead to a narrower bandwidth and thus higher T^, so the 
theory presented here would become sensitive to the hole 
density at higher concentrations. If it is possible to ex- 
tend the control of As antisite concentration pioneered in 
Ref. IT4lto higher Mn densities perhaps still higher Curie 
temperatures may be achieved in very high conductivity 
samples. 

Our studies of the onset of ferromagnetism in a mag- 
netic semiconductor have revealed a mechanism for fer- 
romagnetism in the Mott variable-range-hopping regime 
that is not well described by the commonly-used mean- 
field theory based on the local carrier density. Even 
though carriers are required for formation of the ferro- 
magnetic state, the Curie temperature is almost entirely 
insensitive to the hole density in this regime. The total 
magnetization of the material, however, will be affected 
as the hole density is reduced. Increased disorder leads 
to a larger bandwidth and a lower Curie temperature, 
and eventually to the discontinuous disappearance of the 
ferromagnetic state. Measurements on samples with the 
same experimental Mn concentration and hole concen- 
tration, but very different hopping energies, also indicate 
a strong correlation between hopping energy and Curie 
temperature, and a weak correlation with hole density, 
along the lines of the theory. 
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FIG. 4: (color online) Curie temperature of non-rotated 
Gai-xMuxAs samples with x = 0.0075, 0.01, 0.0125, and 
0.015 plotted as a function of (a) carrier concentration at 300 
K (b) conductivity at 25 K and (c) hopping energy at Tc- 



that the key experimental quantity determining the Curie 
temperature is the hopping energy (or, equivalently, Tq), 
regardless of the microscopic details of the disorder. We 
even find similar results for Ga-rich samples (not shown 
here), which are structurally much less homogeneous [14 
than the As-rich samples described here. As further sup- 
port for this view, samples with W{T = 25K) > 15 meV 
were found to be non- ferromagnetic, whereas all samples 
with W{T = 25^) < 15 meV were ferromagnetic. 

These results currently apply to materials which 
are not at the upper extreme of Curie tempera- 
tures. However, magnetic semiconductors in the low- 
doped regime [25] have longer nonequilibrium carrier and 
spin lifetimes, perhaps making them more suitable for 
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